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RELIABILITY APPROACH TO ROTATING -COMPONENT DESIGN 
by Dimitri B. Kececioglu* and Vincent R. Lalli 
Lewis Research Center 

summary 

The Lewis Research Center, in conjunction with the University of Arizona, has de- 
veloped a design-by-reliability methodology for rotating mechanical components that 
statistically matches the selected material strength to the imposed stresses. Analytical 
and experimental methods have been developed to match materials and loads in this way. 
Design-by-reliability variables are handled as random functions. To illustrate this 
methodology, statistical methods for (1) determining strength distributions for steel ex- 
perimentally, (2) determining a failure theory for stress variations in a rotating shaft 
subjected to reversed bending and steady torque, and (3) relating strength to stress by 
reliability are given in this report. The experimental test machines that provided data 
for the stress used with the design-by-reliability methodology are described. 

A sample rotating component design problem is solved using both a deterministic 
design method and the proposed method. With the proposed method a smaller, lower 
weight shaft can result for a specified life and reliability than that with conventional 
methods. The new design -by-reliability methodology uses the distortion-energy theory 
with statistical fatigue diagrams for optimum shaft design. The conventional determin- 
istic design methodology uses single-valued estimates, which results in unnecessarily 
conservative rotating component design. Other mechanical components, such as springs, 
bearings, gears, clutches, brakes, couplings, linkages, and cams, may also be de- 
signed using the design-by -reliability methodology. 


INTRODUCTION 

Present mechanical design techniques depend heavily on ' ’modifying factors, " such 
as, surface condition, size, temperature, stress concentration, and miscellaneous ef- 
fects, to compensate for our limited knowledge of materials and loads. In the best 

*Professor of Aerospace and Mechanical Engineering, the University of Arizona, 
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cases they are based on empirical data obtained over a long period of time, but in many 
others there are few data, and the selection is based on the designer s feelings at the 
time. A final factor that is often applied and that generally reflects all the uncertainties 
involved is probably influenced as much by the seriousness of the consequences should 
the part fail as by anything else. The resulting design is considerably larger and heav- 
ier than necessary to avoid failure. 

Smaller and lighter designs are important in aerospace technology. Lighter designs 
are possible if the designs are based on the statistical matching of the available strength 
to the imposed stress. Lewis Research Center in conjunction with the University of 
Arizona in 1964 undertook the task of developing a mechanical design methodology of 
this type. The design methodology developed reduced the number of modifying factors 
by providing increased knowledge through laboratory testing methods. It deals with im- 
portant design variables as statistical distributions rather than as single-value esti- 
mates. It also makes provisions to design for a specified life and reliability. To work 
out this methodology the following information is required: 

(1) Statistical strength distributions for selected materials 

(2) Statistical stress variations in a loading configuration 

(3) Statistical methods for relating stress to strength. 

This report presents the analytical, experimental, and application results attained in 

this research. 

A typical rotating component design problem is given in appendix A to illustrate the 
deterministic, single -valued design technique. In appendix B the sample problem is 
worked out using the probabilistic design-by-reliability methodology which shows that 
substantial size and weight savings are accomplished. Very useful probabilistic experi- 
mental data for AISI 4340 steel are also presented. All data referenced in this report 
were taken using U. S. Customary units. To broaden the usefulness of the report, the 
U. S. Customary values were converted to International System(S. I. ) values. The pri- 
mary conversion factors are (1) 1 kilogram is 0.45 pound mass and (2) 1 kilonewton is 
4. 44 pounds force. The S. I. units are the primary reporting values, with the U. S. Cus- 
tomary values given as the second set throughout. 


SYMBOLS 


A, B, C coefficient, staircase test 

Cjj effect of size 

C L effect of type of load 

Co effect of surface finish 
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cr 
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Kb 

K c 

Kf 


K, 


kw 


K s 

Kt 

*v 
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BGR 




BTH 
X B/T 
•^GRTH 


Jfr 


T/B 


a,b,c, 
d, e,f 


specimen radius, cm (in. ) 
specimen major diameter, cm (in.) 
maximum specimen major diameter, cm (in. ) 
critical specimen major diameter, cm (in.) 
diameter 

specimen groove diameter, cm (in. ) 
specimen stress increment, N/m 2 (psi) 

Young’s modulus, N/m 2 (psi) 
gage factor 

moment of inertia, cm (in. ) 
staircase test method stress levels 
polar moment of inertia, cm 4 (in. 4 ) 
theoretical bending concentration factor 
reliability factor 

theoretical fatigue reduction concentration factor 

theoretical keyway concentration factor 

theoretical torsional shear concentration factor 

theoretical part geometry factor 

theoretical constant value 

bending in groove coefficient 

bending to tool holder coefficient 

bending interaction into torque coefficient 

groove to tool holder coefficient 

life cycles coefficient 

tool holder torque coefficient 

torque interaction into bending coefficient 

variation coefficient 

modifying factors in eq. (3) 


breaking tensile load, kN (lbf) 




N a 

N 


N, 


N. 


v 


Ni 


n. 


n 


s 

pdf 

<1 

R c 

R g 

31 

r 

S a 

Sb 

®m 

S n 

S n 




S 5 

s 

s„ 


m 


ultimate tensile load, kN (lbf) 

yield tensile load, kN (lbf) 

bending moment on specimen, m-N (ft-lb) 

number of active arms in strain bridge 

Visicorder divisions 

cycles to failure 

Visicorder 

number of specific cycles of life 
sample size 

th 

number of successes at i level 
number of specimens 
probability 

probability density function 
notch sensitivity factor 
calibration resistance, ohm 
gage resistance, ohm 
reliability 
stress ratio, sVs 

alternating probabilistic Goodman diagram material strength, N/m (psi) 

2 

tensile breaking strength, N/m (psi) 

2 

mean probabilistic Goodman diagram material strength, N/m (psi) 

2 

corrected endurance limit, N/m (psi) 

2 

estimated endurance limit of the rotating beam specimen, N/m (psi) 

2 

tensile ultimate strength, N/m (psi) 

p 

tensile yield strength, N/m (psi) 
stress 

n 

alternating material stress, N/m (psi) 

2 

component strength standard deviation, N/m (psi) 

specific component strength, N/m^ (psi) 

2 

mean material stress, N/m (psi) 
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s n combined stress vector 

s o lowest staircase level material stress, N/m 2 (psi) 

T torque, N-m (lb-ft) 

x, y, z normal random variables 
z standardized normal variable 

atg skewness moment 

<*4 kurtosis moment 

A tolerance, cm (in. ) 

£ difference statistic, N/m 2 (psi) 

0 angle, deg 

M mean estimate 

o standard deviation estimate 

CT x, a alternating bending stress, N/m 2 (psi) 

CT Xj c compression bending stress, N/m 2 (psi) 
c x, m mean bending stress, N/m 2 (psi) 

a x, n notch bending stress, N/m 2 (psi) 

a x, o output bending stress, N/m 2 (psi) 
o Xj t tension bending stress, N/m 2 (psi) 

a alternating radial stress, N/m 2 (psi) 

a y, m naean radial stress, N/m 2 (psi) 

T xz a alternating shear stress, N/m 2 (osi) 

0 

r xz m mean Shear stress, N/m (psi) 

T xz, n notc h shear stress, N/m 2 (psi) 

T xz, o output shear stress, N/m 2 (psi) 

Superscript: 

mean of statistically distributed variables 

BACKGROUND 

The basic probabilistic methodology for designing mechanical components from a 
study of the intersection of their stress and strength distributions was discussed in 



reference 1. Included in that paper was a discussion of Monte Carlo techniques for de- 
termining the stress and strength distributions, given the distributions of the factors 
affecting them. 

Freudenthal (ref. 2) found in his study of structural unreliability that the safety fac- 
tor was a random variable, which was the quotient of strength to stress, where both 
strength and stress were considered random variables. Freudenthal, Garrelts, and 
Shinozuka prepared a comprehensive report (ref. 3) along the same lines which dis- 
cussed in more detail the mathematical techniques required, the appropriate statistical 
distributions involved, and problems that remained. Several example problems in 
structural reliability were worked out, and an extensive bibliography was given. These 
efforts concentrated on fatigue and structural reliability. 

The Battelle Memorial Institute and its Mechanical Reliability Research Center pre- 
sented studies (refs. 4 and 5) that described some of the fundamental problems in me- 
chanical probabilistic design and suggested methods for their solution. 

Mittenbergs (ref. 4) discussed the fundamental aspects of probabilistic design as 
applied to mechanical components. In that paper it was explained that the basic failure 
modes of mechanical components are (1) deformation, (2) fracture, and (3) structural 
instability. The combined effects of the various failure modes must be determined to 
assess the reliability of a mechanical component. The intersection of strength and 
stress distributions was discussed. Reference 5, a summary of 2 years of research, 
contains a thorough discussion of probabilistic design and attempts to quantify the rela- 
tions of various factors on such phenomena as creep and fatigue. An extensive bibliog- 
raphy is included. 

The Illinois Institute of Technology Research Institute conducted an extensive pro- 
gram (ref. 6) in probabilistic design methods. The program was concerned with three 
major areas: (1) The study of prime failure mechanisms in mechanical design (specific 
items included fatigue, wear, creep, and corrosion), (2) the application of failure me- 
chanisms and design information for the reliability assessment of specific mechanical 
components (Parts included were gears, bearings, springs, and shafts.), and (3) the 
determination of mechanical system reliability in terms of individual component 
reliabilities. 

Bratt, Reethof, and Weber (ref. 7) provided a computer approach to the solution of 
the time variant strength distribution case. Haugen (ref. 8) contributed to the under- 
standing of probabilistic loads and their combinations. This work considers two aspects 
of the probabilistic design methodology: (1) the probability theory necessary to perform 
the design function and (2) probabilistic design examples drawn from various engineering 
disciplines. 

Much still remained to be done. Important aspects of the design by reliability 
methodology that remained to be investigated were 
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(1) The development of the design by reliability methodology for specific cases 

(2) The examination of methods to combine functions of random variables as applied 

to design by reliability 

(3) The exploration of current and new methods for determining failure-governing 

stress and strength distributions 

(4) The exploration of current and new methods for determining reliability once the 

failure-governing stress and strength distributions are known 

(5) The development and fabrication of testing machines to provide data for the de- 

sign by reliability methodology 

(6) Carrying out a test program to obtain experimental materials data for use in the 

probabilistic, mechanical, design-by-reliability methodology. 

This research included a literature search, a review of existing theory, design, 
development, and fabrication of test machines ; experimental testing j and the develop- 
ment of computer programs to facilitate the reduction of the experimental data. 

The research effort concentrated on the rotating shaft problem. This methodology 
demonstrates that other mechanical components, such as springs, bearings, gears, 
clutches, brakes, couplings, linkages, and cams can be designed by using reliability 
methods. Six reports (refs. 9 to 15) on this subject describe various aspects of the 
problem of designing rotating components that are subjected to combined -stress fatigue 
by the design-by-reliability methodology. Included are many aspects of estimating dis- 
tributions, estimating stress and strength factors, computing reliabilities by various 
methods, and experimental test procedures. This report describes a probabilistic 
methodology for designing rotating mechanical components using reliability to relate 
stress to strength and the experimental data obtained to support this methodology. 


ANALYSIS 

Determine the Failure Governing Strength 

In the engineering application of the disign-by-reliability methodology, one of the 
most difficult tasks is to determine the actual strength distribution of a specific mate- 
rial in a specific application. Two methods of determining this distribution are syn- 
thesis and direct experimentation. Much work remains to be done in these areas. 

Synthesis . - Synthesis is concerned with finding a statistical function that accurately 
represents the material's strength. Considerable research has been expended by 
Kececioglu (ref. 11) and others to determine which distribution defines the experimental 
data best. The distributions studied are the normal, the lognormal, and Weibull. 
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A procedure to convert readily available tensile strength data, such as that given 
in table A-3 of reference 16 (p. 600) to statistical functions that accurately represent the 
material's endurance strength would be very useful. This table gives tensile yield 
strength S and tensile ultimate strength S^. In many cases where endurance strength 
distribution data for AISI steels are not available, use can be made of the common prac 
tice, based on the work of Lipson and Juvinall (ref. 17, p. 162), to obtain the estimate 
r»f the mean for endurance limit for steels as 

S' = 0. 50 S ut when S ut £ 1378 MN/ra 2 (200 ksi) (1) 

^ = 689. 4 MN/m 2 (100 ksi) when S ut > 1378 MN/m 2 (200 ksi) (2) 

The standard deviation for the normal distribution can be obtained by making use of a 
conclusion based on Kececioglu's and others' work (refs. 9 and 10). These investigators 
have found that a standard deviation of about 7 percent of the unmodified endurance limit 
can be used if little or no test data are available. The final expression 1 of becomes 

(s; ; o.o7 s;). 

Marin's work (ref. 17, p. 127) can be extended to obtain a machined part's distri- 
butional endurance limit S n by considering each parameter in equation (3) as a random 
variable. 


S n * k a k b k c k dWk (3) 

The distribution of products of random variables can be calculated by various methods 
(refs. 18 and 19). Appendix B is an example of such a calculation. In this case (eq. (3)) 
k , k b , and have a value of unity because the experimental results include these 
factors. The reliability factor k c is included by a different method, and is based 
on experimental data rather than on those data obtained from tables as explained pre- 
viously. Appendix C gives the algebra of normal functions method of forming this 
product. 

Another procedure for determining the distribution of the fatigue strength of mate- 
rials is the Monte Carlo technique. References 9 and 20 to 23 disclose a correlation 
between the methods used in this report and these Monte Carlo techniques. Results ob- 
tained from the Monte Carlo computer calculations are within 5 percent of those obtained 
by the algebra of normal functions method when the coefficients of variation are small 
(ref. 9). Considering the experimental errors associated with strength of materials 
measurements, this normal functions method is a simpler and more acceptable proce- 
dure for estimating strength distributions under these conditions. 

formal distribution notation (mean; standard deviation). 
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Djrect experimentation . - In some cases it may not be possible to estimate the dis- 
tribution parameters for the fatigue strength of a material. This almost always is due 
to a lack of statistically meaningful test data or engineering experience with a particular 
material. A laboratory test procedure by which fatigue strength for different materials 
could be obtained is desired. 

The first requirement would be to design suitable test specimens. Specimens must 
be carefully prepared and inspected to insure that they reproduce as many of the essen- 
tial features of the part application as possible (ref. 24). It is often necessary to scale 
the test specimens to retain the same stresses, stress concentration factors, stress 
ratios, and material conditions in a size suitable for testing. A meaningful quantity of 
these test specimens is required. A machine to test the specimens under conditions 
similar to the part application was also needed. 

In this research two types of specimens were used. These are shown in figures 1 
and 2 and are identified as specimen phases I and II. Groups of these specimens were 
tested to failure under fixed simulated environmental conditions. Table I shows the av- 
erage stress levels, stress ratios and number of specimens tested. Figures 3 and 4 show 
the nine distributional stresses as functions of cycles (S-N) diagrams obtained from this 
environmental testing. Carefully prepared specimens were subjected to fixed levels of 
alternating and shear stresses to keep the average stress ratio (alternating/mean) r at 
fixed values. 

Endurance tests were conducted using the "staircase*' method (ref. 26). The re- 
sults are given in table n. This information was used to obtain the distributional Good- 
man diagrams of figure 5. This was accomplished by projecting the endurance strength 
test data to the appropriate stress ratio planes as shown in figure 5(a). Table II shows 
the alternating endurance strength distribution parameters at 25° C from the analysis of 

the test data and the components of the distribution parameters along the various stress 
ratio planes. 

If a probabilistic Goodman diagram for rotational lives shorter than 2. 5xl0 6 cycles 
is desired, diagrams similar to that in figure 5 can be generated by the method shown 
graphically in figure 6. The cycles to failure distributions given in figures 3 and 4 for 
a particular stress ratio and cycles of life are summed. The cumulative histogram is 
formed to show the percentage of specimens failing in each stress cell for a fixed cycle 
life value. From this histogram, the material's combined strength distribution can be 
found (ref. 25). The resulting statistical distribution is the material's strength for a 
fixed number of cycles of life and stress ratio. However, the components of strength 
distribution along the fixed stress ratio plane are needed for specific shaft design prob- 
lems. These strength distribution components can be obtained as follows: 

(1) Obtain values for S a , <r g , r, and Nj from figure 6. 

(2) It is known that 


9 



and 


Sf = ( s a + s m) 


(4) 


1/2 



(5) 


(3) Substitute equation (5) into equation (4) and obtain 



( 6 ) 


(4) Simplify equation (6) and assume that S f and S a are normally distributed and 
that r is a fixed plane. 2 This yields 


and 



(7) 


( 8 ) 


Substituting test values into equations (7) and (8) yields the mean and standard deviation 
components of the strength distribution along r for the probabilistic Goodman diagram. 
It seems that, to establish an acceptable Goodman fatigue strength surface, at least 
four such distributions (r = 00 , 3.5, 0. 83, and 0. 44 at a specified life cycle) are re- 
quired. This, coupled with the static ultimate tensile strength distribution at r = 0, 
gives five sections on the contour to define the probabilistic Goodman strength sui face. 


Consideration of r as a random variable complicates matters considerably. The 
six standard deviations variation should not be greater than ±0.017 rad (±1°) for most 
cases; therefore, r was fixed as a plane. 
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Determine the Failure Governing Stress 


The first problem in determining the failure governing stress is to identify the fail- 
ure governing theory which applies best to the part in service considering its material, 
loads, dimensions, stress concentrations, and the like. Rotating shafts have been ana- 
lyzed using the distortion-energy or von Mises-Hencky failure criterion (refs. 9 and 27). 
This distortion-energy criterion is applicable to ductile steels (ref. 9, pp. 152-154) and 
is theoretically valid for the elastic region only. It is recognized that the steel used in 
this research does not strictly fit the ductile requirement. Although the fracture pattern 
in fatigue of these grooved specimens is ductile; in direct tension the fracture pattern is 
brittle. The elastic region constraint is not strictly met either, as the specimens are 
tested to fracture. The methodology exhibited in this report assumes that the distortion- 
energy theory describes the material stresses that cause failure. Test results indicate 
that other failure theories, combinations of theories, or newly developed theories may 

also accurately reflect the laboratory results (refs. 9, 27, and 28) and they should also 
be studied. 

Figure 7(a) shows the rotating shaft volume element stresses for the general case. 
For the specific case of a rotating shaft subjected to bending and torque stresses only, 
the xz element at the surface simplifies (see fig. 7(b)). Using the stress components 
as defined in figure 7(b) and the von Mises-Hencky failure criterion (ref. 16, p. 154) 
the failure governing alternating stress s a is given by 

s a = a x,a 0) 

and the failure governing mean stress s m by 

s m " T xz, m (10) 

For this stress model, the following design-by-reliability equations can be derived: 



( 11 ) 
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and 


where 



s = = 10. 2 — 
a i/e 



129.7 — 
d 3 


and 


Appendix B gives a detailed example for the use of these equations. 


( 12 ) 


(13) 


(14) 


(15) 


Bridge the Gap by Reliability 


Using reliability to relate material strength to application stresses requiies the 
identification of all failure inodes and the sections where these failures may occur. For 
each such section and failure mode the failure governing strength and stress distribu- 
tions should be determined next, and finally the reliability should be calculated as fol- 
lows: 

A no failure probability exists for the critical failure mode and section when a given 
strength S is not exceeded by a given stress s (ref. 1, p. 555). From figure 8, the 
probability that a stress of value s x exists in the interval s 1 - ds/2 to s x + ds/2 is 
equal to the area of the element ds (A^ as shown in fig. 8) or 

f( Sl )ds = A t (16) 

The probability that S exceeds Sj is equal to the shaded area A 2 in figure 8, or 
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(17) 


^ (S>s l ) = ^°°f(S)dS= A 2 


The probability of no failure is the product of these two probabilities or 


d 31 = f(s^)ds 



f(S)dS 


(18) 


The section’s reliability would then be all probabilities of all strengths being greater 
than all possible values of stress or 



i: 


f(S)dS 


ds 


The section reliability can also be written as (ref. 1, p. 556) 


(19) 


3t = 



dS 


( 20 ) 


Equations (19) or (20) can now be used to calculate the critical section’s reliability 
for any component whose f(s) and f(S) are known. These equations carry limits of inte- 
gration applicable to distributions defined over the interval from to +°°. For func- 
tions defined in different intervals these limits should be replaced by the lowest and 
highest values that can be used. For methods of evaluating equations (19) and (20) see 
references 1 and 25 when neither f(s) and f(S) is normal or lognormal. 

If both f(s) and f(S) are normal, the probability density functions may be expressed 
as 


and 



-l/2[(s-s)/<r ] 2 
e s 


f(S) = 


1 -l/2[(S-S)/crJ 2 


12 1 ) 


( 22 ) 


The critical section reliability is given by the probability that strength is in excess 
of stress or that S - s > 0. Using the designation £ = S - s, the critical section reli- 
ability is given by all of the probabilities that £ > 0. Let f(£) be defined as the differ- 
ence distribution of f(S) and f(s). As f(S) and f(s) are normally distributed, then 
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( 23 ) 


f( 5 ) is also normally distributed and is expressed by (ref. 1, p. 556) 


f(0 = 


x - 1 / 2 [(C -?)/^] 2 

<j^yj2v 


where 


? = S - s 


(24) 


and 


°rM 1/2 (25) 

The critical section reliability would then be given by all probabilities of £ being a 
positive value, hence, 

1 f' e l/2[(C-C)/o ? ] 2 (26) 

a 5 VS J 0 

The relation between C and the standardized normal variate z can be used to 
evaluate equation (26), which is 


z = J^ 


a 


5 


(27) 


The limit of the integrand for £ = 0 is 


z = OjLi = - JL ( 28 ) 

cr^ 

and for £ = 00 is 

z = !L£i = « (29) 

a ? 
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also 


d£ » <7^ dz (30) 

If these conditions are substituted into equation (26), the following result is obtained: 


St = 



-j= e-VaAte 


(31) 


Consequently, the critical section reliability of a component is given by the area under 
the standardized normal probability density function from the value of z = -(£/<r ) to 
z - The value of this area may be obtained from the table of areas under the stan- 
dardized normal density function. 


APPARATUS AND PROCEDURE 

Combining reversed bending with steady torque for fatigue testing of metal speci- 
mens was previously tried by Mabie and Gjesdahl (ref. 29). Their paper explained a 
test machine using the four-square and rotating-beam principles. The four-square 
principle couples torque to a specimen through a connecting shaft, two end gear boxes, 
and a second shaft for torque loading. The rotating beam principle applies reverse 
bending through bearings and a static fixed force. These methods of applying combined 
stresses to metal test specimens without the need to dissipate energy proved satisfac- 
tory. However, this first machine was not able to hold torque or bending stresses con- 
stant, and developed excessive noise and vibration during operation. 

Probabilistic design by reliability applied to space -power system shafts required 
that a test machine capable of providing constant torque and bending stresses for metal 
specimens be developed (ref. 30). Using the proven principles of the Mabie-Gjesdahl 
machine, the University of Arizona designed and built such a machine. The machine 
was carefully reviewed, modified where appropriate, and approved by Lewis. 


Research Machine 


The fatigue machine consists of a two-section, rotating shaft in the front with a test 
specimen locked in the center (fig. 9). The front shaft has a flexible coupling at each 
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end to allow for relatively free deflection when the specimen is loaded in bending. A 
direct 200-watt (7. 5-hp), 1800-rpm motor powers the front shaft. The bending load is 
applied to the specimen by means of two yokes, one on each of two bearings located 
symmetrically about the specimen on two commercial tool holders. Below the shaft 
the yokes are connected by a horizontal link, which concentrates the load at a single 
vertical link in the center. The vertical link is then connected to either a long or a 
short loading lever arm. These loading arms make possible the application of a range 
of bending stresses in the specimen groove by means of pan weights applied at the end 
of the loading arm. The steady torque is applied by means of a commercial Infinit- 
Indexer which is located on the back shaft of the machine which rotates at about 1000 
rpm. Table III summarizes the operational specifications of these research machines. 
The machine is capable of producing, holding, and transmitting to the rotating specimen 
steady torques of up to 605 newton-meters (5400 lb-in. ) and reversed bending moments 
of up to 386 meter-newtons (3450 in. -lb). Specimens with diameters up to 2. 5 centime- 
ters (1 in. ) can be tested in this machine by changing the collets in the specimen holder. 
This machine has the following advantages'. 

(1) A constant bending moment is maintained by directly applied static loads. 

(2) Constant torque is maintained through the use of an Infinit-Indexer. 

(3) Monitoring of specimen strain is reasonably direct. 

(4) Backlash in the four-square loop is eliminated. 

Typical space-power turbine shaft designs were analyzed to make the research data 
directly applicable to aerospace problems. Figures 1 and 2 show the grooved test spec- 
imen details (ref. 31). Bar stock pieces about 15 centimeters (6 in.) long were care- 
fully machined to the dimensions and finish specified in figures 1 and 2. The reverse 
bending moment M, torque T, bending stress cr , and shear stress r xz m are 
identified in figure 1. 


Research Procedure 

Each specimen is installed in the test machine, and the instrumentation is checked 
for zero adjustment and calibration. The appropriate bending moment is applied to the 
specimen by putting weights on the load pan. The torque is applied to the specimen by 
rotating the outer shell of the Infinit-Indexer with a suitable wrench. The timing clock 
is set to zero and the machine is started. When the specimen fractures, a microswitch 
stops the clock and the machine. Figures 1(b) and 10 show the specimen configuration 
and fracture pattern after testing, respectively. 

Strain gages are used on the tool holder to monitor the bending and torque loads. 
These gages are located on the tool holder rather than on the specimen groove. The 
reasons for this are that (1) it is difficult to mount strain gages in the limited space of 
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a. specimen groove and (2) because the specimens are fractured, the gages could also be 
damaged during failure resulting in prohibitive expenses. The positioning and electrical 
circuitry for the strain-gage bridges are given in figure 11. The slip-ring assembly 
is located next to the strain gages. The slip-rings are counterbalanced with an alumi- 
num collar of equal weight and nearly equal dimensions located on the other tool holder 
(see fig. 9). The amplifiers, galvanometers, and recorder are matched units. This 
equipment is used to amplify and record the output from the bending and torque gages. 

Measurement of the nominal bending and shear stresses in the specimen groove are 
required. Because the strain gages cannot be located in the specimen groove, the shear 
and bending stresses in the specimen groove must be determined from the tool holder 
strain gages. To convert tool holder strain data to specimen groove stresses, the fol- 
lowing calibration coefficients were required: 

(1) Actual groove bending stress versus apparent groove bending stress 

(2) Apparent groove bending stress versus apparent tool holder bending stress 

(3) Torque interaction into bending 

(4) Actual groove shear stress versus tool holder shear stress 

(5) Bending interaction into torque. 

Bach machine was carefully calibrated (ref. 12). The relation between the calibra- 
tion variables, in each case, proved to be linear; hence, a slope could be associated 
with each calibration variable. Because for all cases the functional relation began at 
the origin, the slope alone of each curve completely defined the function. These slopes 
have been given calibration coefficient designations and are listed in table IV. The mode 
of use of the calibration results is described in detail in figure 12. Steps 1 and 2 re- 
quire the selection of a stress ratio and a bending stress level. The nominal shear 
stress in the specimen groove to give this ratio is found in step 3. Next, in step 4 the 
nominal bending stress in the groove is converted to output stress in the groove using 

■*BGR and **GRTH' 0n the torc l ue side, step 5 gives the equation for converting the 
nominal shear stress in the groove to shear stress in the tool holder. Also, in step 5 
the shear stress in the tool holder is converted to output stress. Steps 6 and 7 are the 
corrections for interaction between torque and bending. Steps 8 and 9 convert the cor- 
rected output stress to Visicorder divisions, completing the procedure. Each machine 
has the following calibration equation: For bending B 


Nv,B ER ^BGR**GRTH°nom, GR + **T/B T out, Th) ( 3 2) 

g 


17 



and for torque T 


N. 


v, T 


N c N a GR c PBTH . 


ER, 


Jfn 


T out, TH + ^B/T^out, TH 


(33) 


All strain-gage parameters are those of the bending and torque bridges. 

For data reduction a stress ratio and a bending stress level in the groove are se- 
lected, as is shown in blocks 1 and 2 of figure 12. From these two parameters the nom- 
inal shear stress in the groove is determined in block 3. This step uses the distortion- 
energy failure criterion. Next, the nominal bending stress in the groove is converted 
to tool holder stress in block 4. In block 5 the theoretical conversion from the groove 
to tool holder torque and from nominal to output tool holder stress is made for shear. 

In blocks 6 and 7 each tool holder output is corrected for interaction effects. In blocks 
8 and 9 the output tool holder stresses are converted to recorder divisions. Thus, the 
number of recorder divisions necessary to give the required stress ratio and bending 
stress level are determined. 

Next, a sample is run while maintaining the necessary divisions of bending and 
torque as closely as possible. A record of the bending and torsion gage signals is taken 
for every specimen. When the running of the sample lot is complete, the data reduction 
procedure is reversed for each specimen to determine the nominal bending and shear 
stresses and the stress ratio actually achieved. 


EXPERIMENTAL RESULTS 


To reach the primary objectives of this research, it was necessary to obtain dis- 
tributional S-N, endurance, and static tensile strength data. The research consisted of 
running test specimens under selected and controlled conditions for 


r s = =^ = 

S m 


yfi t. 


^ — = 0, 0. 44, 0. 83, 3. 50, *> 


xz, m 


and 


S - (7 

a x, a 


ranging between 172. 4 and 1034. 1 newtons per square meter (25 and 150 ksi). 
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The number of specimens tested for each stress level of an S-N diagram was either 
12 or 18. It was determined that three to six stress levels should be studied to ascer- 
tain the distributional S-N diagram for a particular stress ratio. 

Tables I and II, V to VIII, and figures 3 and 4 present the types of data obtained 
from the research machines and the reduced results. Table I gives stress levels, ra- 
tios, and the cycles-to-failure data (life) for research phases I and n. Table II gives 
the endurance strength distribution results. Table V shows, as an example, how the 
endurance strength data for mean stress ratio of 0. 44 are reduced to obtain their distri- 
bution parameters using the sample data given in figure 13. Table VI gives the grooved 
and ungrooved static tensile strength data as the strength data for r = 0. Tables VII 
and VO summarize the parameters for the normal and the lognormal distributions fit 
to the cycles to failure (life) data. Figures 3 and 4 give the distributional S-N diagram 
presentation of the cycles to failure (life) results given in tables H, VII, and Vin in log- 
normal distribution form. An analysis of the cycles to failure data shows that the dis- 
tributions are log-normally distributed; therefore, they plot as normal distributions on 
log-log scales (ref. 13). Figure 5 gives the probabilistic Goodman diagram for the 
phase I and II data given in table II in normal distribution form. 

Figure 13 is a sample of the data for endurance strength testing. The endurance 
tests were conducted using the ' ’staircase’' method. In the staircase method the bend- 
ing stress is changed in increments of about 20. 7 meganewtons per square meter 
(3. 0 ksi), while maintaining a constant stress ratio. If there is no failure during about a 
2. 5 -million-cycle test, then one increment is added to the bending stress level of the 
next specimen and the test repeated. If a failure occurs, one increment is subtracted 
for the next specimen. The endurance strengths are best represented by normal dis- 
tributions (ref. 10). 


EXPERIMENTAL VERSUS EMPIRICAL S-N DIAGRAM 

If the desired experimental S-N diagram is not available, the designer is forced to 
develop his own based on empirical considerations and long -established rules of thumb. 
To find out how good these rules -of -thumb are, the mean S-N diagrams determined ex- 
perimentally in phase I were superimposed on the conventional thumb -rule developed 
S-N diagram in figure 14. The empirical S-N diagrams were constructed according to 
a procedure recommended by Juvinall (ref. 32, p. 211), knowing the ultimate strength 
of a given steel, as follows; 

For the AISI 4340 steel used in this research, the ultimate strength S ut was 1227. 1 
meganewtons per square meter (178 ksi) for unnotched specimens. The endurance 
strength at 10® cycles of life is given by (ref. 32, pp. 236-256) 
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(34) 


S n C L C D C S 


n i06 


K, 


where 

S n = °- 5 S ut = 613< 6 N / m2 ( 89 ksi ) ( 35 ) 

and where the effect of load type C L equals 1. 0 for bending, the effect of size C D 
equals 0. 9 for diameters between 1 and 5. 08 centimeters (0. 4 and 1 in. ), the effect of 
surface finish C g equals 0. 88 for the fine ground finish, and the fatigue stress concen- 
tration factor Kf is 


% = 1 + qOCj. - l)C g 

where q is the notch sensitivity and Kj. is the theoretical stress concentration factor 
obtained from Peterson's data (ref. 34, p. 49). Hence, for q = 0. 95 and K^. = 1. 42 

Kj = 1 + 0. 95(1. 42)(0. 88) = 1. 351 

Other factors, such as temperature effects, residual stresses, and surface defects, 
are considered to have no effect. Therefore, from equation (34) 

g = (613. 6)(1. 0 )(0. 9)(_0. _ 36Q ( or 52 . 2 ksi) 

n 10 6 1. 351 

The fatigue strength at 10 3 cycles is given by (ref. 32, pp. 260-263) 


where 


n 


= 0. 9 


lO^ 


U ut 

K f 


Kj *= 1 + (K { 


1)Y 


Y= 0.64 


(36) 


and 

S ut = 1227. 1 MN/m 2 (or 178 ksi) 
and where Y is the factor determined by Juvinall. Hence, 
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n l 0 3 


_ J 0^9 )(1227. 1) = 9Q1 7 MN y' m 2 

1 + (0. 351)(0. 64) 


(or 130, 8 ksi) 


In determining the empirical value of the strength at 10 3 cycle life, the unnotched ulti- 
mate strength has been used rather than the notched. To obtain the empirical rule-of- 
thumb based S-N diagram, a straight line is drawn on log- log paper between 901.7 mega- 
newtons per square meter (130.8 ksi) at 10 3 cycles and 360 meganewtons per square me- 
ter (52. 2 ksi) at 10 cycles. This is the empirical line in figure 14, and it represents 
the mean rule-of-thumb cycles-to-failure estimate. The remaining lines shown in fig- 
ure 14 are the best-fit mean cycles-to-failure lines for the four stress ratios researched 
in phase I. 

A comparison of Juvinall's method for constructing an empirical S-N diagram with 
the experimentally determined mean line for a stress ratio of infinity shows the strength 
at 10 cycles is underestimated using Juvinall’s method and the mean value of the ten- 
sile strength of the unnotched specimens. At 10 3 cycles the empirical estimate of the 
mean fatigue strength yielded a value of about 901. 7 meganewtons per square meter 
(130. 8 ksi) compared with 1172 meganewtons per square meter (170 ksi) for a stress 
ratio of infinity, or 75 percent of actual. At the other end of the diagram the empirical 
estimate of the endurance strength yielded a value of 360 meganewtons per square me- 
ter, which is below the experimentally determined mean value of 395 meganewtons per 
square meter for a stress ratio of infinity, or 9 percent lower than actual, thus making 
the empirical estimate conservative. Furthermore, the empirical life of 10 6 cycles is 
high when compared with the experimental value of 700 000 cycles for a stress ratio 
of infinity. Specifically, it is 43 percent higher; consequently, the empirically estab- 
lished endurance life is not conservative. Since Juvinal does not take into account com- 
plex loading involving alternating bending and constant mean torsion, the empirical S-N 
diagram will only be applicable for the stress ratio of infinity. 


DESIGNED EXAMPLE 

A typical rotating component part design problem is given in appendix A. The part 
is the inner shaft of an alternator rotor, shown in figure 15(b). It is designed first using 
a deterministic, single-valued design methodology and modifying factors to reduce the 
material's endurance limit to the case of the problem. A Goodman boundary limit is 
constructed using the modified endurance limit with existing tensile strength data. The 
distortion- energy stress components suitable for this model are theoretically selected. 
Finally, graphical construction methods are used to determine the permissible alter- 
nating and mean stress for this material. The desired shaft diameter is calculated using 
equation (46) of appendix A. The results of this design effort are summarized in 
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table IX, where values are given for the important design parameters. The shaft diam- 
eter following this design procedure came out to be 1. 39 centimeters (35/64 in. ). 

The same shaft is redesigned using the methodology in appendix B. This method 
does not require modifying factors for surface finish, size, or stress concentration as 
the laboratory test program generated data for endurance strengths under combined 
conditions of reversed bending and steady torque for machined and ground test speci- 
mens of about the proper size and with very similar stress concentration. The distri- 
butional Goodman diagram of figure 5(a) is used to determine the failure governing 
strength distribution for this material. The desired shaft diameter is calculated using 
the difference statistics of two normal distributions resulting in equation (63) in appen- 
dix B. The results of this probabilistic design effort are also summarized in table IX. 
The shaft diameter following the probabilistic design procedure came out to be 1. 04 cen- 
timeters (27/64 in. ). This is a 40-percent decrease in the weight of material required 
to fabricate this part. Furthermore, the design is estimated to have a shaft life of 
2. 5x10® cycles and a minimum reliability of 0. 999; that is, it is expected that no more 
than 1 in 1000 such shafts will fail while performing their designed-for mission. 

Using reliability to relate material strength to application stresses requires care- 
ful engineering judgment to assess all failure modes. A given design requires that the 
most likely section to fail and its failure mode be selected and analysis performed at 
that section. In the probabilistic design -by -reliability approach none of the original de- 
sign requirements were reduced. These requirements were met by selecting the most 
probable failure mode and matching it with the material, geometry, and appropriate 
fatigue strength. 

The best designed product is only as good as the people and materials finally used 
to make it. The task of determining that people with the required skills and materials 
of the proper specified quality are used to build the product must be included in this de- 
sign methodology by other means. This task requires participation of many disciplines 
including rigorous quality assurance engineering (ref. 24). There are many important 
design problems that can be solved using the probabilistic design-by-reliability concepts 
developed in this research. Other problems and their solutions using this methodology 
can be found in the literature (refs. 1, 8, 9, 13, 25, and 35 to 37). 


DISCUSSION OF RESULTS 

The research results are summarized in tables I, II, V to VIII and figures 3 and 4 
(distributional S-N diagrams), figure 14 (comparison of the experimental S-N diagram 
with the theoretical), figure 5 and table II (distributional Goodman diagrams), and ta- 
ble IX (comparison of results of conventional design with probabilistic design-by- 
reliability results). 
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The combined- stress fatigue research machines are capable of applying and main- 
taining the desired alternating bending moment and steady torque, based on the low 
achieved ratio of cr r /r. These machines enable the application of about 390 meter- 
newtons (3450 in. -lb) of bending moment and about 600 newton-meters (5400 lb-in. ) of 
torque to the test specimen. 

Probabilistic analysis of the S-N data reveals that the cycles-to-failure distribu- 
tions for various stress ratios are best represented by the lognormal distribution; 
whereas the endurance strength distributions are best represented by the normal dis- 
tribution. 


For phase I results and r = °° the mean cycle life increases substantially as the 
nominal alternating bending stress decreases from about 1170 to 480 meganewtons per 
square meter (170 to 70 ksi). Specifically, the mean of the cycles to failure at 
965 MN/m 2 (140 ksi) is about 2900 cycles, and at 480 MN/m 2 (70 ksi) the mean of the 
cycles to failure is about 187 800 cycles, or a 53 percent increase. The standard devia- 
tion of the cycles to failure at 965 MN/m 2 (140 ksi) is Oj N = 0. 172; at 480 MN/m 2 

(70 ksi) it is cr^ Q g N = 0. 192, or an 11 percent increase on the basis of the logs. It 

must be pointed out that a comparison of the standard deviation on the basis of the 
straight cycles cannot be made because the cycles-to-failure distribution is lognormal 
and the antilog of a^ Q(T ^ does not provide any useful information. The standard devia- 
tion does increase, however, as the alternating stress level decreases although not 
much for r = 00 . 

Figure 3(b), the distributional S-N diagram for r = 3. 5, shows a substantial in- 
crease in the mean cycle-life as the nominal alternating bending stress decreases from 
about 1170 to 480 MN/m 2 (170 to 70 ksi). Specifically, the mean of the cycles to failure 
at 965 MN/m 2 (140 ksi) is about 2100 cycles and at 480 MN/m 2 (70 ksi) is about 98 000 


cycles, or a 50 percent increase on the basis of the logs. The standard deviation of the 
cycles to failure at 965 MN/m 2 (140 ksi) is a logeN = 0. 180; at 480 MN/m 2 (70 ksi) it 
is ^loggN = 286, or a 59 percent increase on the basis of the logs. In this case the 

standard deviation increases substantially as the alternating stress level decreases for 
r = 3. 5 as compared with the situation for r = 


Figure 3(c), the distributional S-N diagram for r = 0. 83, likewise shows a substan- 
tial increase in the mean cycle -life as the nominal alternating bending stress decreases 
from about 830 MN/m 2 (120 ksi) to about 480 MN/m 2 (70 ksi). Specifically, the mean of 
the cycles to failure at 830 MN/m 2 (120 ksi) is about 4300 cycles and at 480 MN/m 2 
(70 ksi) is about 88 200 cycles, or a 36 percent increase on the basis of the logs. The 
standard deviation of the cycles to failure at 830 MN/m 2 (120 ksi) is a log N = 0. 197; at 

480 MN/m 2 (70 ksi) it is a loggN = 0. 197, an increase of only 0. 4 percent. This in- 
crease is even less than that for r = <=°; nonetheless, it is an increase in the standard 
deviation as the alternating stress level decreases. 
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Figure 3(d), the distributional S-N diagram for r = 0. 44, shows a smaller increase 
in the mean cycle life as the nominal alternating bending stress decreases than for the 
other stress ratios, but the range of the alternating bending stress is also smaller, 
namely, from about 550 MN/m 2 (80 ksi) to about 415 MN/m 2 (60 ksi). Specifically, the 
mean of the cycles to failure at 550 MN/m 2 (80 ksi) is about 19 100 cycles; at 
415 MN/m 2 (60 ksi) it is about 130 500 cycles, or an increase of 20 percent. The stand- 
ard deviation of the cycles to failure at 550 MN/m 2 (80 ksi) is <J l0 g^ N = 0- 245; at 
415 MN/m 2 (60 ksi) it is a lQg N = 0. 329, or an increase of 34 percent. Thus, the 

standard deviation likewise increases as the alternating stress level decreases. 

Conclusions similar to these may be drawn from a study of figure 4, for phase II 
results. 

A study of figure 3 shows that the following maximum nominal alternating bending 
stress levels are permissible for phase I results: For r = «>, max = 1700 MN/m 
(170 ksi), for r = 3. 5, S a max = 1170 MN/m 2 (170 ksi), for r = 0. 83, S a> max = 

860 MN/m 2 (125 ksi), andfor r = 0. 44, S a max = 535 MN/m (78 ksi) to avoid impos- 
ing bending and shear stresses in the vicinity of the yield strength of the material. 

A study of figure 4 shows that the following maximum nominal alternating bending 
stress levels are permissible for phase II results: For r = °°, S a ^ max - 895 MN/m 
(130 ksi); for r = 1. 06, S a max = 830 MN/m 2 (120 ksi); for r = 0. 40, max = 

375 MN/m 2 (55 ksi); for r’= 0. 25, max = 3 10 MN/m (45 ksi); andfor r = 0. 15, 

S = 225 MN/m 2 (37 ksi). 

5 A comparison of figure 3 and a study of table III, all for phase I results, indicate 
that, for the same alternating stress level, as the stress ratio decreases from ®° to 
0. 44 the mean life also decreases. In particular, at a nominal alternating stress level 
of 480 MN/m 2 (70 ksi) the following mean lives are exhibited: For r = 00 about 187 800 
cycles, for r = 3. 5 about 98 100 cycles, for r = 0.83 about 88 200 cycles, andfor 
r = 0 . 44 about 46 600 cycles. It is apparent that there is a significant reduction in the 
mean life as the shear stress, which is superimposed onto the bending stress, is in- 
creased, that is, as the stress ratio is decreased. From r = to 3. 5 the mean life 
(in cycles) is reduced by 48 percent; from r ~ 3. 5 to 0. 83 the mean life is reduced by 
10 percent; and from r = 0. 83 to 0. 44 the mean life is reduced by another 47 percent. 

A comparison of figure 4 and a study of table H for phase II indicate that similar 

conclusions may be drawn for phase n results. 

A further comparison of figure 3 and a study of table II indicate that, as the stress 
ratio decreases from <*> to 0. 44, the standard deviation of the cycles to failure in- 
creases in general. In particular at a nominal alternating stress level of 480 MN/m 
(70 ksi) the following standard deviations are exhibited: For r = », cr log N = 0. 192; 
for r = 3. 5, a lQg N = 0. 286; for r = 0. 83, a j Qg N = 0. 197; and for r = 0. 44, o log N = 
0. 274. These data show that standard deviation increases in terms of log 10 cycles to 
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failure, as the stress ratio decreases. A similar trend may also be observed for other 
alternating stress levels. 

A comparison of figures 3(a) and 4(a) shows that at an alternating bending stress 
level of 480 MN/m^ (70 ksi) and r = 00 , the mean cycles to failure decreases from about 
187 800 cycles for phase I to about 12 800 cycles for phase II results, or a 93 percent 
decrease. This sharp reduction in the mean life results when the groove geometry is 
so. altered that its theoretical stress concentration factor increases from 1. 42 to 2. 34. 
Consequently, great care must be exercised not to introduce such high stress concen- 
trations. 

A comparison of figures 3 and 4 indicates that, as the theoretical stress concentra- 
tion factor is increased from 1. 42 to 2. 34, the standard deviation of the cycles -to- 
failure decreases because of the greater concentration of stress in a smaller volume of 
groove material and the higher actual stresses in the groove, 

A study of figure 14 indicates that for r = 00 Juvinall ! s method for constructing an 
empirical S-N diagram underestimates the fatigue strength at 10 3 cycles of life by about 
25 percent because the experimental alternating bending strength at 10 3 cycles is 1170 
MN/m 2 (170 ksi) instead of the 901. 7 MN/m 2 (130. 8 ksi) predicted by Juvinall's em- 
pirical method. It appears that either one or both factors that determine the empirical 
fatigue strength at 10 3 cycles (eq. (36)) ought to be modified. Specifically, if it is as- 
sumed that the fatigue stress concentration factor K[ is the best available to date, a 
correction for the 0. 9 factor would be in order such that the new factor for this mate- 
rial and geometry would be 1170/1176 (170/178) = 0. 955. 

Figure 14 indicates that the mean endurance strength at 10® cycles of life is also 
underestimated by Juvinall’s empirical method. Specifically, the experimental value is 
395 MN/m 2 (57. 3 ksi) as compared with 360 MN/m 2 (52. 2 ksi) for the empirical esti- 
mate. A study of equations (34) and (35) indicates that the potential contributors to this 
difference might be primarily the factor of 0. 5 in equation (35) and the factor C D or 
size effect in equation (34). It may be assumed that, because of lack of good recent data 
on C D , the correction may best be made on the factor 0. 5, which could be modified to 
395/360 (57 300/52 200) = 0. 549. Therefore, equation (35) would now be 


S n = °- 549 S ut 


(37) 


It is of interest to compare the cycles of life at the knee of the empirical S-N dia- 
gram with that of the experimental S-N diagram for r = It may be seen that the ex- 
perimental endurance life is 700 000 cycles as compared with 10® cycles for the theoret- 
ical of Juvinall, or 30 percent lower. Furthermore, even greater reductions in the 
endurance life result from the superimposing of a constant torsional shear stress onto 
the alternating bending stress; namely, the life at the knee is 300 000 cycles for r=3. 5, 
270 000 cycles for r = 0. 83, and 550 000 cycles for r = 0. 44. 
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A study of table 11(a) indicates that, as the stress ratio decreases from 00 to 0. 45, 
the endurance strength under combined bending -torsion decreases from a nominal al- 
ternating bending stress of 395 MN/m 2 (57. 3 ksi) to 341. 9 MN/m^ (49. 6 ksi), or a 
13 percent decrease, which points out the significant effect on the mean endurance 
strength of the superimposition of increasing torsional shear stress onto the alternating 
bending stress at 2. 5xl0 6 cycles of life. 

A study of table 11(b) indicates that, as the stress ratio decreases from 00 to 0. 15, 
the endurance strength under combined bending -torsion decreases from a nominal al- 
ternating bending stress of 232. 3 MN/m 2 (33. 7 ksi) to 157. 9 MN/m^ (22. 9 ksi), or a 
34 percent decrease. This decrease, again, points out the significant effect on the mean 
endurance strength of the super imposition of increasing constant torsional shear stress 

onto the alternating bending stress at 2. 5xl0 6 cycles of life. 

The distributional Goodman diagrams of figure 5, for 2. 5X10 cycles of life, show 
that there is little effect on the alternating bending strength mean at stress ratios 
greater than 1. 0, due to the superimposition of constant torsion onto alternating bend- 
ing stress. However, there is an increasingly sharp reduction in the alternating bend- 
ing strength as the stress ratio decreases below 1. 0. It appears that the added shear 
stress is asserting its influence when its nominal stress level is higher than about 
230 MN/m 2 (33 ksi). Furthermore, figure 5(a) indicates that there is a substantial in- 
crease in the standard deviation of the failure governing strength distribution as the 
stress ratio decreases below 1. 0 and that the standard deviation reaches a maximum 
at a stress ratio of about 0. 30. Thereafter, it drops sharply and converges to its min- 
imum at r = 0, which is the standard deviation of the static ultimate tensile strength 

distribution. 


CONCLUDING REMARKS 

The abscissas in figure 5 show the static ultimate tensile strength of the ungrooved 
specimens because this is the most commonly available static strength information. A 
Goodman diagram analysis using the static ultimate strength distribution of the grooved 
specimens has been programmed. Further, nominal alternating bending strength has 
been used throughout to represent the alternating strength. There are two reasons for 
this: (1) The actual bending and shear stresses at the base of the test specimen groove 
are not known because the exact fatigue stress concentration factor for this material 
under combined bending torsion is not known. (2) Normally in design, the nominal 
bending and shear stresses are calculated first, and, to determine the corresponding 
fatigue strength distribution, the designer need only enter figures 5(a) or (b) with these 
nominal stresses. 
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The design application presented in this study provides a vivid illustration of the 
great value of this research. It shows how to calculate the reliability of shafts subjected 
to a variety of loading conditions by the design-by-reliability methodology. The calcu- 
lated reliability provides an explicit quantitative measure of the design integrity that no 
other method can provide. The 40-percent weight savings achieved by the design-by- 
reliability methodology over that of the conventional design methodology in the design 
example for the same reliability dramatizes the value of the methodology and of the ex- 
perimental, distributional S-N, and Goodman diagram design data generated by this 
research. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 1, 1974, 

120-60. 
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APPENDIX A 


DETERMINISTIC SHAFT DESIGN 


The shaft shown in figure 15 supports an alternator rotor . It is to be made of 
AISI 4340 steel, heat treated to a Brinell hardness of 323 to 370 (Rq 35 to 40) with a 
ground finish. The rotor will subject the shaft to a radial bending load of 607. 4 newtons 
(136. 8 lb) and a torque of 112 newton-meters (1000 lb-in.) when delivering full-load 
electrical power. The spline drive may exert a maximum misalinement radial force of 
13. 3 newtons (3 lb). The temperature to be experienced by the shaft will be ^294 to 300 K 
(70° to 80° F). The shaft will rotate at 12 000 rpm and have a life of 2. 5XlO b cycles 
with a reliability of 0. 999. 

What should the shaft’s diameter be using the conventional deterministic design 


methodology? 

Solution: Various methods can be used for machine design problems in which the 
shaft diameter is unknown. For aerospace applications weight, life, and reliability are 
important considerations. Shigley(ref. 16, p. 482 ) has shown that under these conditions 
it is wise to use the von Mises-Hencky-Goodman method. This method requires knowl- 
edge of material strengths and application stresses. Table A-3 of reference 32 (p. 564) 
gives the tensile properties for 2. 54-centimeter (1-in. ) round, drawn condition, 

277 Brinell hardness material as Sy = 1103 MN/m 2 (or 160 ksi) and S ut = 1289 MN/m 

(or 187 ksi). 

According to Juvinall (ref. 32, pp. 236-256) the endurance limit S n may be found 
using modifying factors as follows: 


n 10 6 


S n C L C D C S 

Kc 


(34) 


where 

= 0, 50 S ut = 644. 6 MN/ra 2 (or 93. 5 ksi) 

when S t = 1379 N/m 2 (200 ksi; ref. 32, p. 211), and where the type of load C L is 1.0 
for reverse bending (ref. 32, p. 231), the effect of size is 0. 9 for specimens up 
to 5. 08 centimeters (2 in. ; ref. 32, p. 232), the effect of surface C g is 0. 88 for a 
ground finish (ref. 32, p. 234), and the fatigue stress concentration factor K f is 
(ref. 32, p. 257) 


Kj = 1 + (Kj. - l)qCg 


(38) 
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The maximum theoretical shear and bending stress concentration factor at this shoulder 
are (ref. 34, pp. 75-76) 


K g = 1. 19 and = 1. 30 

for a shaft around 1. 27 centimeters (1/2 in. ) in diameter with a 0. 318-centimeter 
shoulder radius. What to do when stress concentration is caused by a combined bending- 
torsion load has not been adequately investigated, but it is reasonable that the combined 
effect would be no greater than the product of each factor; therefore, 

K t = (1. 19)(1. 30) = 1.55 

Notch sensitivity q, based on data (ref. 33, pp. 296-298) for a notch radius of 0. 318 
centimeter (0. 125 in.) and S yt = 1289 MN/m 2 (187 ksi), is 0. 90. Therefore, the maxi- 
mum fatigue stress concentration factor from equation (38) is 

K f = 1 + (1. 55 - 1)(0. 90) (0. 88) = 1. 43 


From equation (34) 


S = 644. 6(1. 0)(0. 90 ) 0 X 88) = 357 ! ^^2 j or gL 8 ksi ) 
n 1.43 


To satisfy the reliability specification, an additional modifying factor K c (ref. 16, 
p. 169) must be introduced. From Shigley (ref. 16, p. 169) K c is determined as 

K c = 1 - 0. 008 D = 0.75 


The resulting endurance limit now becomes 


S n, corr = K c S n = (0 - 75 H 357 - *) = 268. 9 MN/m 2 (or 39. 0 ksi) 

The bending moment M at the shoulder is 14. 1 meter-newtons (126 in. -lb). The 
bending stress a at the shoulder is calculated from 2 

" A, d 


Me _ 14. 1 _ 144 

1 jrd 3 /32 d 3 



(39) 


Sparentheticals contain U. S. Customary values. 
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The torque T is constant over the shaft length at 113 newton-meters (1000 lb-in. ). The 

torsional stress t v „ „ in the design section of the shaft is 
xz. m 


T xz, m 


Tc 

J 


113 _ 576 

7rd 3 /16 d 3 



(40) 


From the von Mises-Hencky failure criterion, given by equations (5), (9), and (10), 
the failure governing mean shear stress -r is maximum at the outer fiber of the 

AZ ? III 

shaft at which point the contribution due to the vertical shear of the specimen is zero 
and therefore, is not included in the analysis. Thus, the von Mises-Hencky stress 
components using equations (9) and (10) become 


i -a =!« 
a x, a o 

d 6 



(41) 


and 


m 


= J~3r =998 
xz. m „ 

d d 


for 


8830) 


(42) 


and equation (5) becomes 


s a 

r = — - = 0. 144 
s m 


(43) 


A conventional Goodman diagram can be prepared from these data. The alternating 
and mean stresses are plotted as ordinate and abscissa, respectively, in figure 16. De- 
fining stress-to-failure as material strength makes it possible to use the Goodman 
boundary to represent the strength of this material as modified by environmental factors. 
The top portion of the Goodman boundary (line 1) is obtained by drawing a straight line 
between the ordinate at = 269 MN/m 2 (39 ksi) and the abscissa at S u{; = 1289 MN/m 2 
(187 ksi). The right boundary (line 2) is defined by a straight line joining the ordinate 
at S v = 1103 MN/m 2 (160 ksi) with the abscissa also at S y = 1103 MN/m 2 (160 ksi). A 
third straight line (line 3) is constructed to define S„ and S , for this case, with 
tan 9 = 0 . 144 x 6 = tan 40. 8. (The factor of 6 takes into account the fact that the abscissa 
scale is 6 times that of the ordinate scale. ) The intersection of lines 1 and 3 define the 
point S a = 109. 6 MN/m 2 (15. 9 ksi) and S m = 762. 5 MN/m 2 (110. 6 ksi). 
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In conventional design a margin of safety of at least 100 percent is used to insure a 
no failure design. This requires that s be constrained as follows: 

cl 


_ . a _ 109. 6 

a_ 2 2 


(“ 1 r) 


(44) 


Substitution of equation (41) into equation (44) yields 


144 

d 3 


< 54.8 


(or -^2 < 7. 95X10 3 ] 


(45) 


Then 


d 3 > Jii- = 2.627 
54.8 


for 1280 = 0.161^ 
\ 7950 J 


(46) 


or 


d > 1. 38 cm (or 0. 544 in. ) 


Based on these results a shaft diameter of at least 1. 39 centimeters (35/64 in. ) 
should be recommended. The inputs and results of the conventional deterministic design 
methodology are summarized in table EX for purposes of comparison with the results of 
the design-by-reliability methodology, which is discussed in appendix B. 
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APPENDIX B 


SHAFT DESIGN BY RELIABILITY 


Now let the alternator shaft of appendix A be designed using the reliability approach 
to rotating component design. 

Solution: To design by reliability, data in the form of the distributional Goodman 
diagram given in figure 5(a) must be used. The design -by-reliability methodology has 
been detailed in a previous paper (ref. 1). The shaft is to be made of AISI 4340 steel, 
heat treated to Rockwell C hardness 35 to 40 or Brinell hardness 323 to 370. Figure 5(a) 

provides the fatigue strength data required. 

Since the research has generated data for the strength of the material desired under 
combined conditions of reversed bending and steady torque for various stress ratios, 
the modifying factors required by Juvinall in appendix A, namely, C D and C g are al- 
ready included in the data. The reliability factor K c is introduced differently in the 
design-by-reliability methodology, as will be shown later. The difference in K f for 
this shaft and the research specimens is negligible. All other modifying factors are 
taken to be unity as in appendix A. 

The bending moment at the shoulder is taken to be normally distributed 1 with pa- 
rameters (14. 1; 1. 58) meter-newtons ((126; 11. 4) in. -lb). The failure governing stress 
s at the shoulder is calculated from equation (13). The distribution parameters of s , 

el 41 

taken to be normal, are found using equation (13) and the method of binary synthesis of 
normal distributions (ref. 38, pp. 279-282 and table 9), such that 


s a =10.2 


M 

d 3 


(47) 


and 


where 



(48) 


^■Normal distribution notation (mean; standard deviation). 


32 



= .3 d 2 <r d = 3 


d 2 (0. 015 d) = 0.045 d 3 


and where a d is taken to be 1/6 of the total specified tolerance, or 0. 015 d, M = 14. 1 
meter-newtons (126 in. -lb) and a M = 1. 58 meter-newtons (11. 4 in. -lb). (This has to be 
specified. ) Substitution of these quantities into equations (47) and (48) yields 


10.2 


14. 1 
d 3 




(49) 


and 



10. 2 J (14. 1) 2 [0. 045 d 3 1 2 + d 6 (1. 58) 2 
d 6 + (0.045 d 3 ) 2 


1/2 


d 3 




(50) 


The distribution parameters of s m , taken to be normal, are found similarly from equa- 
tion (42) as follows: 


s 


m 



(51) 


and 


where 



(52) 


T = 113 N-m (or 1000 lb-in. ) 
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and 


ct t = 9 N-m (or 80 lb-in. ) 

Substitution of these quantities into equations (51) and (52) yields 


« = 8 . 83 = MS or 8830^ 

d 3 d 3 [ d 3 


( 53 ) 


and 


2 

8.83 l (113) 2 fo. 045 d ^1 + d 6 (9) 2 
d 6 + (0.045 d 3 ) 2 


S - d 3 


1/2 


^ far ) (54) 

d 3 \ d 3 ' 


From equations (50), (54), and (43) 


144 


r = 


a d 


3 3 


= 0. 144 


S m 998 
d 3 


(55) 


(It is assumed that the variability of r will not affect the reliability significantly. ) The 
two stresses, s a and s m , must now be synthesized into the failure governing stress 
distribution f(s f ), along the stress ratio line r = 0. 144, equations (7) and (8), as fol- 
lows: 


s f 



(56) 


and 



( 57 ) 
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In equation (56) s a is obtained from equation (49) and r = 0. 144. And in equation (57) 

°s, a obtained from equation (50). Substitution of these in equations (56) and (57) 
yields 


il/2 


s f = 144 (1 + 

3 3 


d 3 \ (0. 144) 2 i 


= 144 — (7. 016) = 1010 — i- 


d 3 


d 3 


or 


9016 
d 3 


(58) 


Similarly, 


a = 14.7 JL[i + 1 

f d 3 \ (0. 144) 2 


a/2 


= 14. 7 — (7. 016) - 103 — 


d 3 


d 3 


or 


910 
d 3 ' 


(59) 


Thus the mean and the standard deviation of the failure governing stress distribution 
have been determined. 

The next step is to determine the failure governing strength distribution parame- 
ters. Since the design strength data presented in figure 5(a) apply to this case for 
r = 0. 144. Figure 5(a) yields 


S f = 1213 MN/m 2 (176 ksi) 


and 


= 46. 4 MN/m 2 


(6. 73 ksi) 


If the failure governing strength distribution f(S f ) and the failure governing stress dis- 
tribution f(sj) are taken to be normal, then these two distributions can be coupled to ob- 
tain with the reliability of these shafts as follows (ref. 1, pp. 556-557): 


31 


r+oc 

C/cr 


K/« c ) « z,dz 


( 60 ) 
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where 



( 61 ) 


As the specified shaft reliability is 0. 999, the value of z(f) may be found to be -3. 09 
from the standardized normal distribution probability tables. Consequently, using this 
value and the values of the previously determined parameters of f(s f ) and f(S f ), equa- 
tion (61) yields 


-3.09 


1010 


- 1213 


d 3 


103 


d 3 


I 2 ? 

+ (46. 4 r 


1 1/2 


(62) 


Rearrangement of equation (62) yields 

d 6 - 1.6885 d 3 + 0.6332 = 0 (63) 

and solving for the roots of equation (63) yields a mean shaft diameter of d = 1. 04 
centimeter (0. 411 in. ). Consequently, the use of d = 1. 07 centimeter (27/64 in. ) would 
be recommended, as compared with a diameter of d = 1. 39 centimeters (or 35/64 in. ) 
required by the conventional, deterministic design. 

The inputs, outputs, and results for both design methodologies are given in ta- 
ble IX. An a priori reliability is designed into the shaft through the design-by- 
reliability methodology, which results in a better than 40-percent saving in material for 
this case. Furthermore, the design integrity is known in terms of a reliability of 0. 999, 
such that on the average no more than one in 1000 such shafts will fail while performing 
a mission of 2. 5xl0 6 cycles duration. If a higher reliability is desired, the correspond- 
ing value for - (S/o^) can be found, and the required shaft diameter recalculated by the 
design-by- reliability methodology. Since all design variables are explicitedly brought 
into this methodology and incorporated into the failure governing stress and strength dis- 
tributions to enable the determination of the designed-in reliability, the design can be 
optimized from the point of view of shaft geometry, loads, and material strength 
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characteristics. This example thus vividly illustrates the great value of the probabil- 
istic design-by-reliability methodology of designing rotating shafts subjected to com- 
bined reversed bending and steady torque loads. 
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APPENDIX C 


THE ALGEBRA OF NORMAL FUNCTIONS 

The algebra of normal functions is a useful method of performing mathematical 
operations with statistical functions (ref. 8). This method is primarily intended for use 
in the probabilistic design -by -reliability approach. All variables are assumed to be 
normally distributed; the mathematical operations for the sum, difference, product, 
and quotient of two normal random variables are also assumed to be normally distri- 
buted. (This is only valid for the sum and difference of two normally distributed random 
variables. ) Under these conditions, an algebra can be developed to facilitate mathemat- 
ical operations required in design. The mathematical operations of this algebra are 
summarized in table X. Some variables that are important in the design-by-reliability 
method are not normally distributed. To handle these design variables other methods, 
such as the Monte Carlo (ref. 39 or 40) and the system moments (ref. 40, pp. 228-235), 
are available. 
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TABLE L. - COMBINED BENDING AND SHEAR STRESS LEVELS 
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m « jvZi ? in 


and cr 


s„ 


= Vs» T 


Not applicable. 


xz, m 
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TABLE II. - SUMMARY OF ENDURANCE STRENGTH RESULTS FOR 
DISTRIBUTIONAL GOODMAN DIAGRAM 


[2. 5x10® cycles of life, AISI 4030 steel, R c 35 to 40; grooved specimens.] 


Stress 

ratio, 

r 

Nominal alternating bending 
strength 

Nominal shear strength 3 

Nominal combined strength 
vector b 

T xz, 

m 

Standard de- 
viation, 

s a 

Standard de- 
viation, 

a 


Standard de- 
viation, <j c 

MN/m 2 

ksi 

MN/m 2 

ksi 

xz, 

ai 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

Phase I 




00 

395. 0 

57.3 

20.0 

2. 9 

0 

0 

0 

0 

395.0 

57.3 

20.0 

2.9 

3.5 

379. 9 

55. 1 

25. 5 

3.7 

108.2 

15.7 

7.6 

1. 1 

395.0 

57. 3 

26.9 

3.9 

1.0 

384. 7 

55.8 

22, 8 

3. 3 

391.6 

56.8 

22.8 

3.3 

543. 9 

78.9 

32.4 

4.7 

. 45 

341.9 

49.6 

25. 5 

3.7 

759. 7 

110.2 

56. 5 

8, 2 

833. 5 

120. 9 

62.0 

9.0 


Phase II 

00 

232.3 

33. 7 

18. 6 

2. 7 

0 

0 

0 

0 

232.3 

33.7 

18.6 

2. 7 

1.06 

220.6 

32.0 

19. 3 

2. 8 

208.2 

30.2 

17.9 

2.6 

303. 3 

44. 0 

26.2 

3.8 

.40 

209.6 

30.4 

20.0 

2. 9 

523.9 

76.0 

49.6 

7.2 

564.6 

81. 9 

53.8 

7. 8 

.25 

166. 1 

24. 1 

11.7 

1.7 

665.3 

96. 5 

47.6 

6.9 

686. 0 

99. 5 

48.9 

7. 1 

. 15 

157.9 

22.9 

4. 8 

.7 

1052. 7 

152.7 

32.4 

4.7 

1064. 4 

154. 4 

32.4 

4,7 


T = a 

and 

<7 -_J_ a 

xz, m ^5 - 

— r- — ‘ 

xz, m y 3 r a a 

_ / A 

k l/2 

/ 

\ = S a 1 + — 

\ rV 

1 

and <j q = cr 0 1 1 4 - — 

s v M p 
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TABLE III. - OPERATIONAL SPECIFICATIONS FOR THE COMBINED- STRESS 


FATIGUE RESEARCH MACHINES 


Component 

Description 

Test machine a 

Accommodate a specimen rotating at 1800 rpm; produce and hold a steady torque 
and a reversed bending moment; holding chuck, 2. 5-cm-diam maximum; simple 
design employing "off-the-shelf 1 ’ components. 

Loading mechanism for 
steady torque 

Simple device to produce, hold, and transmit desired steady torque of 605 N-m 
(5400 lb-in.) to test specimen. 

Loading mechanism for 
reversed bending 

Simple device to produce a reversed bending moment of 386 m-N (3450 in. -lb) 
while specimen is rotating 

Test specimen 

AISI 4340 steel, condition C-4; MIL-S-5000B, certification of chemical and phys- 
ical properties; uniform quality, same heat and processing, heat treat to Rockwell 
C 35 to 40 as per MIL-H-6875 with minimum tempering temperature of 540° C; 
inspection as per MIL-I-6868; D = 1. 87 cm, d = 1. 27 cm, r = 0. 38 cm; 1^=1. 45; 
K =122 with 2. 5 by 0. 32 cm keyway. 

s — 

Instrumentation - 
Strain gages 
Channels 
Slip rings 
Amplifier 
Recorder 

To obtain dynamic and static strain measurements in bending and torsion. 

To handle at least six sets of strain gage output simultaneously. 

To transfer strain gage data to amplifier while specimen is rotating. 

To amplify static and dynamic output from strain gages. 

To produce a permanent record of amplified strain gage outputs. j 


| \ - " 

a Three such machines have been built and are being used to obtain distributional combined stress fatigue 


data for the steel specimens specified above. 


TABLE IV. - MODE 5 a CALIBRATION COEFFICIENTS 


Machine 

a b 

Calibration coefficient , 

Machine 

speed, 

rpm 

Bending in 
groove, 

■*BGR 

Groove to 
tool holder, 

•*GRTH 

Tool holder 
torque., 

Jf T 

Torque inter- 
action into 
bending, 

•*t/b 

Bending inter- 
action into 
torque, 

•*B/T 

1 

1.3081 

0.0149 

0.8753 

-0. 0462 

0.0477 

1779 

2 

1. 3081 

* 0145 

.8950 

.0463 

.0555 

1775 

3 

1.3081 

.0170 

.8748 

.0542 

. 0645 

1778 


^ode 5 operation refers to the computer program code that identifies the proper set of 
calibration coefficients to be used with each set of machine data obtained during a speci- 


fic calendar period. 

k These calibration constants are traceable to the National Bureau of Standards standards. 
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TABLE V. - TYPICAL E DURANCE STRENGTH 


DATA PROCESSING 21 


(AISI 4340 steel; Rockwell C, 35 to 40; stress ratio 
ratio, 0. 44; grooved specimen. See fig. 1. 

Test was considered a success if specimens did 
not break in 2. 5x10** cycles (see fig. 21).] 


Alternating 

i 

n i 

in i 


strength, 

(level) 

(successes) 


S a 






MN/m 2 

ksi 





370 

53.7 

3 

3 

9 

21 

350 

50. 8 

2 

8 

16 

32 

330 

47.9 

1 

7 

7 

7 

310 

45.0 

0 

3 

0 

0 

Totals: 


n = 21, 

A = 32, 

B = 66 


a A calculation of the mean and standard deviation 


using the staircase method is 
dj = 20. 0 MN/m 2 (2. 9 ksi) 

S Q = 310 MN/m 2 (45. 0 ksi) 


^S a = s 0 + d i(| + ^ = 3 1. 0 + 2. 0^ + 350 MN/m 2 


(50. 8 ksi) 


= 1.620 d| 
b a 


21 X 66 - (32) 2 

(2D 2 


1 ^ E 2 A2 + Q .029j = 
o. 029J = 


= 1. 620 x 2. 0 


= 29 MN/m 2 (4. 2 ksi) 


A summary of the mean and standard deviation en- 
durance strength estimates are given in table II. 




►f* 

05 


TABLE VII. - LIFE DISTRIBUTION PARAMETERS FOR PHASE I RESULTS 


Stress 

ratio, 

r 

Alternating 
stress, s 

Sample 

size, 

n 

Life, kilo- 
cycles 

Skewness, 

“3 

Kurtosis, 

a 4 

K-S test 
difference, 

Life, 

kiiocy< 

L°g e 

:les 

Skewness, 

«3 

Kurtosis, 

“4 

K-S test 
difference, 
D 

(a) 

MN/m 2 

ksi 

N f 


i) 

(a) 

Nf 

“N f 



Normal dist 

ributions 

Log 

-normal dis 

tributions 


aO 

696. 3 

101.0 

12 

2.8 

0.53 

-1.344 

3. 561 

0. 175 

7.907 

0.228 

-1.586 

4. 258 

0. 207 


551.5 

80.0 

18 

9.0 

1.02 

-.247 

1. 995 

.099 

9. 102 

.114 

-.407 

2. 123 

.092 


475. 7 

69.0 



22.2 

3. 82 

-.042 

1. 855 

,082 

9.992 

. 176 

-.265 

1.945 

.094 


396.4 

57.5 



78.0 

12.55 

.873 

3.060 

, 184 

11. 252 

. 154 

.575 

2.697 

. 159 


355.0 

51. 5 



162.0 

34.39 

-.616 

2.201 

. 174 

11, 971 

.235 

-.906 

2.707 

. 198 

3.5 

730. 8 

106.0 

12 

1. 5 

0.28 

-0.222 

1.725 

0, 113 

7. 262 

0. 198 

-0. 435 

2.006 

0. 128 


550.0 

80. 5 

18 

6.2 

1.01 

,797 

3.397 

. Ill 

8. 721 

. 157 

.403 

2, 869 

.086 


403.3 

58. 5 

18 

39.6 

12.99 

1. 876 

7.241 

.239 

10. 545 

.287 

.686 

4,407 

. 184 


358. 5 

52.0 

18 

74.2 

19.53 

. 147 

2.057 

. 132 

11. 180 

.274 

-.342 

2. 501 

.084 

0. 83 

537. 7 

78.0 

12 

6.6 

1.01 

9. 893 

2. 493 

b 0. 296 

8. 778 

0. 146 

0. 742 

2. 295 

c 0. 280 


444.7 

64. 5 

18 

20.5 

5.59 

.239 

2.257 

. 119 

9. 890 

.281 

-.171 

1.916 

. 120 


365. 4 

53.0 

18 

61.0 

11. 13 

-.428 

2. 182 

, 109 

11.001 

.195 

-.703 

2.457 

. Ill 


313.7 

45.5 

18 

127.6 

21.23 

. 171 

1.942 

. 141 

11.743 

.167 

-.030 

1. 769 

.147 

0. 44 

334.4 

48. 5 

18 

53.6 

15.47 

1.667 

4.833 

0.272 

10. 857 

0.238 

1.210 

3. 899 

0.214 


310.2 

45.0 

18 

83.7 

26.35 

-.315 

2. 408 

.067 

11. 277 

.374 

-1.091 

3.747 

.115 


289. 5 

42.0 

18 

142.6 

41. 59 

,393 

3.057 

. 148 

11. 825 

.305 

-.457 

3. 182 

. 118 


a The Kolmogrov-Smirnov test rejection or acceptance is based on the absolute value comparison Of maximum D with critical D 
given as follows: 


Level of 

Critical D 

significance, 

value 

Qt 

n - 12 

n = 18 

0.01 

0.338 

0.278 

.05 

.375 

1 

.309 


b Maximum D for normal distributions. 
c Maximum D for log-normal distributions. 





TABLE VIIL - LIFE DISTRIBUTION PARAMETERS FOR PHASE H RESULTS 


[Sample size, 35.] 


Stress 

ratio, 

r 

Alternating 
stress, S 

Life, kil 

Locycles 

Skew- 

ness, 

"3 

Kurto- 

sis, 

a 4 

K-S test 
difference, 
D 
(a) 

Chi -squared test 

Life, log e 
kilocycles 

Skew- 

ness, 

a 3 

Kurto- 

sis, 

«4 

K-S test 
difference, 
D 
(a) 

Chi -squared test 

Nf 

*N f 

Degrees 
of free- 
dom 

Value 

(b) 

Degrees 
of free- 
dom 

Value 

(b) 

MN/m 2 

ksi 

Nf 

^f 


Normal distributions 





Log-normal distributions 



oo 

750. 8 

108. 9 

3 300 

360 

-0. 163 

2.317 

0. 094 

1 

0. 248 

8. 095 

0. 112 

-0. 374 

2. 408 

0. 096 

1 

2.844 


634.9 

92. 1 

6 400 

610 

-.325 

1.672 

. 134 

2 

6. 151 

8. 760 

.097 

-. 402 

1.684 

. 137 

2 

8. 534 


507.4 

73,6 

14 200 

1 370 

447 

2.654 

. 108 

2 

. 880 

9.554 

. 100 

664 

2.936 

. 128 

1 

.255 

1. 06 

728.7 

105.6 

4 500 

650 

0.252 

2. 786 

0. 109 

1 

1.265 

8. 407 

0. 145 

-0. 141 

2,991 

0.087 

1 

0.288 


587. 4 

85. 2 

10 800 

1 210 

.200 

2. 187 

. Ill 

2 

.726 

9.281 

. 112 

.016 

2.089 

. Ill 

1 

.728 


447.4 

64.9 

24 300 

2 760 

.757 

2.095 

. 139 

1 

1. 095 

10.093 

. 110 

.513 

2.751 

. 117 

2 

1. 595 


277. 1 
1 

40.2 

224 000 

83 280 

.650 

3.143 

. 112 

1 

2.453 

12.252 

.380 

. 175 

2.347 

. 114 

2 

8. 188 

0.40 

342.6 

49. 7 

59 000 

9 610 

0. 969 

3.626 

0. 154 

1 

3.029 

10.973 

0. 155 

0. 567 

3. 126 

0. 122 

0 

(c) 


275, 8 

40.0 

210 700 

36 740 

.735 

3.364 

. 103 

1 

.631 

12.244 

. 169 

.313 

2.722 

.077 

2 

0. 892 


242,0 

35. 1 

332 700 

110 670 

1. 478 

5.138 

c .201 

1 

2.928 

12.670 

.296 

.682 

3.237 

. 142 

2 

1.206 

0.25 

' 273.0 

39.6 

57 120 

14 560 

0. 383 

2. 220 

0. 140 

1 

2.606 

10.921 

0.257 

-0.075 

2.352 

0. 094 

2 

2. 862 


| 241.3 

35.0 

95 700 

22 490 

. 597 

i 2.790 

. 099 

2 

3. 106 

11. 442 

.231 

. 135 

2.291 

.068 

2 

.528 

0. 15 

221. 3 

32. 1 

148 000 

31 800 

0. 526 

2.758 

0. 093 

1 

1.017 

11. 883 

0.212 

0. 101 

2.300 

0. 095 

2 

4.410 


188.9 

27.4 

272 100 

64 500 

. 578 

3.881 

. 106 

1 

3, 875 

12.487 

.240 

-.237 

3.278 

.086 

1 

3.011 


179.2 

26.0 

504 100 

156 480 

1. 560 

6. 194 

. 194 

0 

<d> 

13,091 

.277 

. 722 

3.272 

e .147 

1 

1. 226 


a The critical diameter D fir is 0. 224 for a 0. 05 level of significance and a sample size of 35. 

The accept or reject critical value for chi-squared test at 0,05 level of significance is 3. 841 for 1 degree of freedom and 5.991 for two degrees of 


freedom. 

c Maximum D for normal distributions. 

j 

a Not applicable. 

e Maximum D for log-normal distributions. 
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TABLE IX. - ALTERNATOR ROTOR SHAFT DESIGN SUMMARY OF DESIGN 


BY RELIABILITY AND CONVENTIONAL DESIGN METHODOLOGIES 


Parameter 

Conventional design ap- 
proach values 

Design -by -reliability approach values 

C L 

C D 

KC 

C S 

K f 

S n 

corr 

T 

M 

s a 

s m 

S m 

d 

& 

1.0 

0.90 

0. 75 
0.88 

1. 43 

644. 6 MN/m 2 (93. 5 ksi) 
268. 9 MN/m 2 (39.0 ksi) 
113 N-m (1000 lb-in.) 
14, 1 m-N (126 in. -lb) 

i* MN/m 2 f 1 * 28 tail 
d 3 _d 3 

998 MN/m 2 r 8 - 83 ksil 

d 3 [_d 3 

109.6 MN/m 2 (15.9 ksi) 
762.5 MN/ra 2 (110.6 ksi) 
1. 39 cm (35/64 in. ) 

<c) 

Not required 
Not required 
Not used 
Not required 
Not required 
Not required 
Not required 

a (112; 9) N-m (1000; 80 lb-in.) 

(14.1; 1.58) m-N (126; 11.4 in. -lb) 

f 14 i; ^Wm 2 !” 1 - 28 ; 13 ksi 

^d 3 d 3 _ d 3 d 3 J 

l 98 ; 91, 5 1 MN/m 2 8 ' °' 81 ksi 

[_d 3 d 3 , )_d 3 d 3 j 

b ( 12 13.3; 46.4) MN/m 2 (176.0; 6.73 ksi) 

1.07 cm (27/64 in.) 

0. 999 


a Normal distribution notation (mean; standard deviation). 
b These are the parameters of the failure governing strength distribution 
f(Sj) along r = 0. 144. 

c The reliability obtained by the approach of Shigley (ref. 16, p. 169) is not 
theoretically applicable to this complex case; hence, the reliability is 
really not known for the conventional methodology. 




TABLE X. - ALGEBRA OF NORMAL FUNCTIONS SUMMARY 

[The expressions for mean and standard deviation for sums of differences of 
normal variates are exact. The remainder of the expressions are approx- 
imate only, and they are good approximations when the coefficient of varia- 
tion V, of x and y are small, i. e. , V = ip/n) < 0. 10 and 

v y = («y p y ) < 0. 10.] 
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(a) Dimensions, finish, forces, and stresses. 



(hi Before and after tests. 

Figure 1. - Phase I specimen geometry. 


1.21 -cm diam.-\ x Keyway, 0.475 





■7. 62 cm 

—15. 24 cm 



1.88 diam. 



Figure 2. - Phase II research specimen geometry. 




Nominal alternating stress, ksi 


1.265 cm i 0.498 In.) dram* 



(a) Stress ratio, op. 

20x10^ 



<0 Stress ratio, 0.83. 

20X10 2 



Log io cycles to failure 


(d) Stress ratio, 0.44. 

Figure 3. - Cycles-to-fallure distributions of A IS 1 4340 steel; Rockwell C 35 to 40 phase I qrooved 
specimens. 3 
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Endurance strength, ksi 



Nominal alternating stress, ksi 



— | 50 

— 40 

— 30 

— 20 


— 40 

— 30 

— 2D 


40 
H30 
— *20 


— 30 

— 26 
— 22 
— 18 


Stresses above this level are 

h not feasible for design f Cycles-to-failure 

I distribution 

Mean-"' 

Endurance strength distribution 
y i I 1 1 1 ( 1 1 l i i 1 1 i l 1 1 I i ( i M il I i I 1 1 1 111 


: 




I0 4 I0 5 

Cycles to failure 

i i i i L 


to 6 10 7 

J_J 1 


3.0 3.4 3.8 4.2 4.6 5.0 5.4 5.8 6.2 6.6 7.0 

log^Q cycles to failure 

(e) Stress ratio, 0. 15, 

Figure 4. - Cycles-to-failure distributions. A IS I 4340 steel, R c 35 to 40 steel grooved specimens. 
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Endurance strength, ksi 



Nominal alternating strength, ksi 


r s 0. 45 


<a) Phase I results,- theoretical part geometry factor, 1. 42. (See table III and fig, 1 . ) 




0 100 200 300 400 500 600 700 800 900 100U 1100 1200 

Nominal mean stress, MNjm 2 


60 80 100 120 140 160 180 

Nominal mean stress, ksi 


(b) Phase II results; theoretical part geometry factor, 2. 34, (See table IV and fig, 2 . ) 

Figure 5. - Distributional Goodman strength diagram. 2. 5xlG 6 Cycles of life; AISI 4340 steel, Rockwell C hardness, 33 to 40. 


r Cumulative histogram 

/ r Cumulative distribution 

I \ 

! \ | r Strength distributions 

r ~ ~ \ i / at 10 5 cycles of life 


- Cycles -to-fai lure 
distributions 


v Mean value for fixed r. 


Cycles to failure 

Figure 6. - Cumulative stress ~to-fal1u re distribution to strength distribution 
conversion technique. 












Specimen holder 


Infinit-indexer 


C$-*f»0778 


Slip rings 


Loading linkage for bending 

■ ..,11 


Specimen 


Brushes 


Figure 9. - Closeup of combined-stress - fatigue research machine. 





(a) Stress ratio, 0.83; alternating bending stress, 770 MN/m 
(111.7 ksi); mean shear stress, 534 MN/m 2 (77.5 ksi). 


(b) Stress ratio, alternating bending stress, 790 MN/m 2 
(114.6 ksi); mean shear stress, 0. 


Figure 10. - Closeup of fractured section at notch. 
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To sliprings 



vClamping device Strain gage schematic 

\ (rotating) 



Figure 11. - Test setup for bending calibration. 
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Figure 12. - Calibration flow chart showing use of calibration results given in table IV to 
assure proper specimen environment. 











Nomina) alternating 
stress, ksf 


• Successful test, i.e. , specimen 
did not break In 2.5x10 6 cycles 



Number of specimens 

Figure 13. - Endurance strength data. Stress ratio, 0. 44; 5AE 4340 
steel; Rockwell C hardness, 35 to 40; grooved specimens. 


20xIQ 2 



10 3 lO 4 IO 5 10 6 10 7 


Cycles to failure 

Figure 14. - Experimental and empirical S-N diagrams based on phase 1 results. 
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P ■ 607. 4 N 
<136. a lb) 


Spline 

centerline 


Pj = 13. 3 N 
<3 lb! 

f 2. 218 cm 
1.874 in. I 



w- 1.13 kg 
12.5 lb) 

1 


1.868 cm 
(0.735 in.) 


R L -46. 

(10.51b) 


tIO. 318-cm ~ 
L 6N 1 ni?-in.l rad, <K8-inJ rad/ 


24. 76 cm 
(9.75 in. ) 


i_ 23,5 cm 
rT" (9. 25 in. 


-Rotor key centerline 





R 


1 

u 


/ 


-Bearing centerline 

£)} I ■ 112 N-m 
u (1000 lb-in. > 


-R D = 585,2 N (131.81b) 


2,54 cm (1 in. 1 


(a) Inner shaft details. 



(bi Bending moment diagram. 

Figure 15. - Alternator inner rotor shaft 



iiii i l i 1 I 

0 24 48 72 96 120 144 168 192 

Mean stress or strength, ksi 

Figure 16. - Conventional combined -stress fatigue diagram show- 
ing modified Goodman line, Soderberg line, and stress ratio 
line for deterministic shaft design example. 
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NASA- Langley, 1975 


E-5528 



